package ‘openxlsx’ was built under R version 3.6.2
Merge templates with 14C, C, and N data
Radiocarbon analyses for the 2001 samples were not run originally, but were completed on archived samples in 2020.
Fig. 16. Litter incubation \(\Delta\)14C-CO2 (2019)
Caption: Mean \(\Delta\)14C-CO2 for each site. Error bars show min and max of duplicate incubation samples. a) Data shown by site, without litter depth, b) Data shown by depth of litter layer, binned by climate zone.
Fig. 16. Mean profile \(\Delta\)14C for 2001 samples
Caption: Mean \(\Delta\)14C by depth for each site in 2001. Error bars show ±1 standard deviation, solid vertical line shows \(\Delta\)14C of the atmosphere in the year of sampling.
Fig. 16. Profile \(\Delta\)14C for 2009 samples
Caption: Profile \(\Delta\)14C by depth for each site in 2009. Solid vertical line shows \(\Delta\)14C of the atmosphere in the year of sampling. Error bars not shown as only a single replicate profile was analyzed per site.
Fig. 16. Mean profile \(\Delta\)14C for 2019 samples
Caption: Mean \(\Delta\)14C by depth for each site in 2019. Error bars show ±1 standard deviation, solid vertical line shows \(\Delta\)14C of the atmosphere in the year of sampling.
Fig. 16. Mean profile \(\Delta\)14C for 2001 and 2019 samples
Caption: Mean \(\Delta\)14C by depth for each site in 2001 and 2019. Error bars show ±1 standard deviation. Vertical lines show \(\Delta\)14C of the atmosphere in 2001 (solid) and 2019 (dashed).
Fig. 16. \(\Delta\)14C-CO2 of 2019 bulk soil incubations
Caption: \(\Delta\)14CO2 by depth for each site in 2019. One rep from GRrf 10-20 (the 10-20 cm increment sample from the cold granite site) is strongly depleted relative to the other rep: \(\Delta\)14C-CO2 = -396.7, -23.5. The highly depleted sample has been excluded for display reasons.
Fig. 16. \(\Delta\)14C-CO2 of 2001 bulk soil incubations
Caption: \(\Delta\)14CO2 by depth for each site in 2001. Note that some sites only have two depth increments. Similar to the 2019 dataset, one of the GRrf reps from the deepest depth increment was strongly depleted: \(\Delta\)14C-CO2 = -469.1. Both points have been excluded for display reasons.
Fig. 16. \(\Delta\)14C-CO2 of 2001 and 2019 bulk soil incubations
Caption: \(\Delta\)14CO2 by depth for each site in 2001 and 2019. Different depth increments were sampled in 2001 and 2019. Points are the mean of laboratory duplicates; error bars are the measured values of each duplicate. Granite/cold point exlcuded for display reasons as it is strongly depleted.
Fig. 16. Respiration data from incubations of 2019 and 2001 bulk soils.
Caption: Points show measured CO2 production of laboratory duplicates as cumulative fluxes or daily flux rates by depth, lines show the means, and the ribbon represents the range.
Fig. 16. \(\Delta\)14C of 2019 bulk soil incubations and corresponding bulk soil
Caption: \(\Delta\)14C of bulk soil and respired CO2 by depth for each site in 2019. Error bars show one standard deviation for bulk soil, points show mean of three replicate profiles for bulk soils and single observations for respired CO2.
Fig. 16. \(\Delta\)14C of 2001 bulk soil incubations and corresponding bulk soil
Caption: \(\Delta\)14C of bulk soil and respired CO2 by depth for each site in 2001. Points show mean of three replicate profiles for bulk soils and mean of laboratory duplicates for respired CO2. The incubated soil samples are a composite made by homogenizing subsamples from each of the three replicate profile samples by depth. Error bars show one standard deviation for bulk soil and the measured values from laboratory duplicates of the incubated composite samples.
$emtrends
PM d14c_mean.bulk.trend SE df lower.CL upper.CL
AN 0.514 0.144 42 0.224 0.804
BS 0.873 0.169 42 0.532 1.214
GR 0.952 0.136 42 0.679 1.226
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
AN - BS -0.3593 0.222 42 -1.621 0.2481
AN - GR -0.4387 0.198 42 -2.221 0.0792
BS - GR -0.0794 0.216 42 -0.367 0.9287
P value adjustment: tukey method for comparing a family of 3 estimates
$emtrends
eco d14c_mean.bulk.trend SE df lower.CL upper.CL
pp 0.779 0.145 42 0.487 1.071
rf 0.777 0.317 42 0.137 1.417
wf 0.606 0.147 42 0.308 0.903
Confidence level used: 0.95
$contrasts
contrast estimate SE df t.ratio p.value
pp - rf 0.00222 0.349 42 0.006 1.0000
pp - wf 0.17373 0.206 42 0.841 0.6797
rf - wf 0.17150 0.350 42 0.491 0.8762
P value adjustment: tukey method for comparing a family of 3 estimates
Fig. 16. Regression of 2019 bulk soil incubations and corresponding bulk soil \(\Delta\)14C
Caption: Regressions of \(\Delta\)14C of bulk soil and respired CO2 by depth for each site in 2019. Error bars show one standard deviation for bulk soil, points show mean of three replicate profiles for bulk soils and single observations for respired CO2.
Fig. 16. Time series of \(\Delta\)14C by depth, as measured
Caption: Points show mean of three profile replicates for 2001, 2009, and 2019 samples. Error bars show ± 1 standard deviation of the mean (only a single profile was analyzed in 2009). Stars show litter incubation \(\Delta\)14C-CO2 for 2019 samples as a point of reference.
Soils collected in both the 2001 and 2009 sampling campaigns were sampled by horizon, but the depth intervals differed between the two sampling years. In 2009, full profiles were excavated for each site, as opposed to the shorter profiles collected in 2001 from the GR and AN sites. Radiocarbon was measured on all three replicate profiles at each site for the 2001 samples, but only for one of the replicate profiles at each site in 2009, e.g. ANpp rep2, etc.
In order to compare the radiocarbon profiles between 2001, 2009, and 2019 we first interpolated both radiocarbon and carbon stock data at 1 cm intervals for each site in the datasets from each year. The carbon-stock-weighted radiocarbon values for any given target depth interval can then be calculated as a simple sum of the product of the carbon weight of each 1 cm increment (relative to the total carbon stock of the target depth interval) and its radiocarbon value. A monotonic cubic spline fit with Hyman filtering was used for the carbon stock interpolation (Wendt and Hauser 2013), and a mass-preserving spline was used to fit the radiocarbon data (Bishop, T.F.A., McBratney, A.B., Laslett, G.M., (1999) Modelling soil attribute depth functions with equal-area quadratic smoothing splines. Geoderma, 91(1-2): 27-45).
Fig. 16. Time series of bulk soil \(\Delta\)14C by 2001 depths (2001, 2009, 2019 samples)
Caption: Points for 2001 samples show the mean \(\Delta\)14C values at the measured depths. Points for 2009 and 2019 samples are spline-fitted estimates of \(\Delta\)14C predicted for the same depth intervals as measured in 2001. Error bars show ± 1 standard deviation of the mean of three replicate profiles for 2001 and 2019 samples (only a single profile was analyzed in 2009).
Fig. 16. Time series of bulk soil \(\Delta\)14C by depth (splined to 2019 depths)
Caption: Points for 2019 samples show the mean \(\Delta\)14C values at the measured depths. Points for 2001 and 2009 samples are spline-fitted estimates of \(\Delta\)14C predicted for the same depth intervals as measured in 2019. Error bars show ± 1 standard deviation of the mean of three replicate profiles for 2001 and 2019 samples (only a single profile was analyzed in 2009). NB: Only two depth intervals were measured at the cool and cold andesite sites (max depth of 27 and 28 cm, respectively), so linear extrapolation (using the slope of the last 1cm spline-fitted depth increment) was used to extend the profiles to 30 cm.
$`10`
$`20`
$`30`
Fig. 16. Change in \(\Delta\)14C of bulk soil (panel a) and respired CO2 (panel b) over time relative to the atmosphere
Caption: Points for 2019 samples show the mean \(\Delta\)14C values at the measured depths. Points for 2001 and 2009 (bulk only) samples are spline-fitted estimates of \(\Delta\)14C predicted for the same depth intervals as measured in 2019. Error bars for bulk samples in panel (a) show ± 1 standard deviation of the mean of three replicate profiles for 2001 and 2019 samples (only a single profile was analyzed in 2009); error bars for incubation samples in panel (b) show the values of the two reps, while the point represents the mean. NB: Only two depth intervals were measured at the cool and cold andesite sites (max depth of 27 and 28 cm, respectively), so linear extrapolation (using the slope of the last 1cm spline-fitted depth increment) was used to extend the profiles to 30 cm.
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
Fitting mass preserving splines per profile...
|
| | 0%
|
|========================================================================================| 100%
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2001 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| 3.776 | -57.040 | 64.592 | 0.987 | BS-AN |
| 43.034 | -17.782 | 103.850 | 0.201 | GR-AN |
| 39.258 | -19.742 | 98.258 | 0.239 | GR-BS |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2001 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| 47.757 | 4.712 | 90.803 | 0.028 | BS-AN |
| 45.292 | 2.247 | 88.338 | 0.038 | GR-AN |
| -2.465 | -44.225 | 39.295 | 0.988 | GR-BS |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2001 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| 55.395 | 17.734 | 93.055 | 0.003 | BS-AN |
| 66.204 | 28.544 | 103.865 | 0.001 | GR-AN |
| 10.810 | -25.726 | 47.346 | 0.742 | GR-BS |
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2001 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -71.522 | -121.112 | -21.931 | 0.004 | rf-pp |
| -69.488 | -117.597 | -21.379 | 0.004 | wf-pp |
| 2.033 | -47.557 | 51.623 | 0.994 | wf-rf |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2001 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -30.492 | -77.815 | 16.831 | 0.260 | rf-pp |
| -33.125 | -79.035 | 12.785 | 0.190 | wf-pp |
| -2.633 | -49.956 | 44.690 | 0.989 | wf-rf |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2001 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -19.910 | -71.507 | 31.686 | 0.605 | rf-pp |
| -9.482 | -59.538 | 40.574 | 0.884 | wf-pp |
| 10.429 | -41.168 | 62.025 | 0.869 | wf-rf |
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2001 inc PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -8.178 | -63.772 | 47.417 | 0.923 | BS-AN |
| -3.095 | -58.689 | 52.500 | 0.989 | GR-AN |
| 5.083 | -50.512 | 60.678 | 0.969 | GR-BS |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2001 inc PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| 23.755 | -39.544 | 87.055 | 0.600 | BS-AN |
| -44.282 | -110.671 | 22.108 | 0.224 | GR-AN |
| -68.037 | -134.426 | -1.648 | 0.044 | GR-BS |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2001 inc PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| 41.442 | -56.220 | 139.105 | 0.523 | BS-AN |
| -61.708 | -164.138 | 40.721 | 0.288 | GR-AN |
| -103.151 | -205.580 | -0.721 | 0.048 | GR-BS |
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2001 inc ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -54.988 | -96.644 | -13.332 | 0.010 | rf-pp |
| -34.104 | -75.760 | 7.552 | 0.118 | wf-pp |
| 20.883 | -20.773 | 62.539 | 0.416 | wf-rf |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2001 inc ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -21.298 | -100.627 | 58.032 | 0.766 | rf-pp |
| 6.974 | -68.664 | 82.612 | 0.968 | wf-pp |
| 28.272 | -51.058 | 107.601 | 0.629 | wf-rf |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2001 inc ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| 8.875 | -106.631 | 124.380 | 0.978 | rf-pp |
| 61.977 | -48.153 | 172.107 | 0.333 | wf-pp |
| 53.103 | -62.403 | 168.608 | 0.471 | wf-rf |
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2019 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -41.756 | -76.898 | -6.613 | 0.018 | AN-GR |
| -9.522 | -44.665 | 25.621 | 0.779 | BS-GR |
| 32.233 | -2.909 | 67.376 | 0.077 | BS-AN |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2019 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -54.667 | -106.932 | -2.401 | 0.039 | AN-GR |
| -16.478 | -68.743 | 35.788 | 0.714 | BS-GR |
| 38.189 | -14.077 | 90.455 | 0.183 | BS-AN |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2019 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -35.544 | -84.884 | 13.796 | 0.191 | AN-GR |
| -25.789 | -75.129 | 23.551 | 0.406 | BS-GR |
| 9.756 | -39.584 | 59.096 | 0.875 | BS-AN |
$`40`
| 30-40 cm | ||||
|---|---|---|---|---|
| 2019 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -31.478 | -89.509 | 26.553 | 0.380 | AN-GR |
| -19.733 | -77.765 | 38.298 | 0.677 | BS-GR |
| 11.744 | -46.287 | 69.776 | 0.869 | BS-AN |
$`50`
| 40-50 cm | ||||
|---|---|---|---|---|
| 2019 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -46.722 | -107.064 | 13.620 | 0.151 | AN-GR |
| -27.589 | -87.931 | 32.753 | 0.498 | BS-GR |
| 19.133 | -41.209 | 79.475 | 0.712 | BS-AN |
$`60`
| 50-60 cm | ||||
|---|---|---|---|---|
| 2019 bulk PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -9.289 | -72.136 | 53.558 | 0.928 | AN-GR |
| -7.044 | -69.891 | 55.802 | 0.958 | BS-GR |
| 2.244 | -60.602 | 65.091 | 0.996 | BS-AN |
$`70`
NULL
$`80`
NULL
$`90`
NULL
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2019 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -31.500 | -65.613 | 2.613 | 0.074 | wf-rf |
| 14.222 | -19.890 | 48.335 | 0.559 | pp-rf |
| 45.722 | 11.610 | 79.835 | 0.007 | pp-wf |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2019 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -35.556 | -82.679 | 11.568 | 0.165 | wf-rf |
| 35.944 | -11.179 | 83.068 | 0.159 | pp-rf |
| 71.500 | 24.376 | 118.624 | 0.002 | pp-wf |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2019 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -58.167 | -94.986 | -21.348 | 0.002 | wf-rf |
| 10.767 | -26.052 | 47.586 | 0.748 | pp-rf |
| 68.933 | 32.114 | 105.752 | 0.000 | pp-wf |
$`40`
| 30-40 cm | ||||
|---|---|---|---|---|
| 2019 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -67.222 | -107.865 | -26.579 | 0.001 | wf-rf |
| 14.578 | -26.065 | 55.221 | 0.648 | pp-rf |
| 81.800 | 41.157 | 122.443 | 0.000 | pp-wf |
$`50`
| 40-50 cm | ||||
|---|---|---|---|---|
| 2019 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -59.322 | -105.522 | -13.122 | 0.010 | wf-rf |
| 28.344 | -17.855 | 74.544 | 0.294 | pp-rf |
| 87.667 | 41.467 | 133.866 | 0.000 | pp-wf |
$`60`
| 50-60 cm | ||||
|---|---|---|---|---|
| 2019 bulk ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -59.622 | -114.836 | -4.409 | 0.033 | wf-rf |
| -27.478 | -82.691 | 27.736 | 0.440 | pp-rf |
| 32.144 | -23.069 | 87.358 | 0.330 | pp-wf |
$`70`
NULL
$`80`
NULL
$`90`
NULL
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2019 inc PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -3.917 | -40.999 | 33.165 | 0.959 | BS-AN |
| 7.583 | -29.499 | 44.665 | 0.857 | GR-AN |
| 11.500 | -25.582 | 48.582 | 0.705 | GR-BS |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2019 inc PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -16.133 | -76.600 | 44.333 | 0.768 | BS-AN |
| -24.507 | -87.925 | 38.911 | 0.582 | GR-AN |
| -8.373 | -71.791 | 55.045 | 0.937 | GR-BS |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2019 inc PM | ||||
| PM.diff | PM.lwr | PM.upr | PM.p.adj | Pairs |
| -36.657 | -86.659 | 13.346 | 0.170 | BS-AN |
| -21.283 | -68.959 | 26.392 | 0.490 | GR-AN |
| 15.373 | -34.629 | 65.376 | 0.706 | GR-BS |
$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| 2019 inc ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -19.217 | -54.416 | 15.983 | 0.357 | rf-pp |
| -2.367 | -37.566 | 32.833 | 0.983 | wf-pp |
| 16.850 | -18.350 | 52.050 | 0.447 | wf-rf |
$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| 2019 inc ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -32.333 | -87.289 | 22.623 | 0.303 | rf-pp |
| -48.617 | -101.015 | 3.782 | 0.071 | wf-pp |
| -16.283 | -71.239 | 38.673 | 0.724 | wf-rf |
$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| 2019 inc ECO | ||||
| ECO.diff | ECO.lwr | ECO.upr | ECO.p.adj | Pairs |
| -14.457 | -68.364 | 39.450 | 0.766 | rf-pp |
| -21.933 | -73.332 | 29.465 | 0.520 | wf-pp |
| -7.477 | -61.384 | 46.430 | 0.930 | wf-rf |
$ANpp
$ANpp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| ANpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -108.646 | -161.123 | -56.17 | 0.005 | 2019-2001 |
$ANpp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| ANpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -39.132 | -121.572 | 43.307 | 0.258 | 2019-2001 |
$ANpp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| ANpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 19.62 | -60.937 | 100.178 | 0.536 | 2019-2001 |
$ANrf
$ANrf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| ANrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -39.662 | -97.2 | 17.876 | 0.116 | 2019-2001 |
$ANrf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| ANrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -23.451 | -62.372 | 15.471 | 0.151 | 2019-2001 |
$ANrf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| ANrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 3.602 | -44.344 | 51.548 | 0.826 | 2019-2001 |
$ANwf
$ANwf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| ANwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 2.096 | -54.912 | 59.104 | 0.924 | 2019-2001 |
$ANwf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| ANwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 7.388 | -22.778 | 37.553 | 0.534 | 2019-2001 |
$ANwf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| ANwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 5.913 | -13.918 | 25.744 | 0.454 | 2019-2001 |
$BSpp
$BSpp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BSpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -34.257 | -85.03 | 16.516 | 0.134 | 2019-2001 |
$BSpp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BSpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -4.4 | -43.862 | 35.063 | 0.772 | 2019-2001 |
$BSpp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BSpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -16.345 | -82.688 | 49.998 | 0.532 | 2019-2001 |
$BSrf
$BSrf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BSrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 14.282 | -37.804 | 66.367 | 0.489 | 2019-2001 |
$BSrf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BSrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -8.656 | -67.505 | 50.193 | 0.704 | 2019-2001 |
$BSrf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BSrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 27.438 | -46.279 | 101.155 | 0.36 | 2019-2001 |
$BSwf
$BSwf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BSwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -36.239 | -83.527 | 11.049 | 0.1 | 2019-2001 |
$BSwf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BSwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -64.075 | -115.738 | -12.412 | 0.026 | 2019-2001 |
$BSwf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BSwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -113.138 | -153.195 | -73.08 | 0.001 | 2019-2001 |
$GRpp
$GRpp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GRpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -51.541 | -120.524 | 17.443 | 0.107 | 2019-2001 |
$GRpp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GRpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 36.232 | -32.197 | 104.661 | 0.215 | 2019-2001 |
$GRpp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GRpp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 2.096 | -49.559 | 53.75 | 0.916 | 2019-2001 |
$GRrf
$GRrf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GRrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -4.878 | -56.92 | 47.164 | 0.808 | 2019-2001 |
$GRrf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GRrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 9.203 | -13.93 | 32.335 | 0.331 | 2019-2001 |
$GRrf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GRrf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 6.556 | -18.916 | 32.027 | 0.514 | 2019-2001 |
$GRwf
$GRwf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GRwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -89.004 | -114.543 | -63.464 | 0.001 | 2019-2001 |
$GRwf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GRwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -65.737 | -141.304 | 9.829 | 0.073 | 2019-2001 |
$GRwf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GRwf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -65.759 | -140.676 | 9.159 | 0.071 | 2019-2001 |
$AN
$AN$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| AN 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -47.196 | -96.931 | 2.54 | 0.061 | 2019-2001 |
$AN$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| AN 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -16.142 | -50.009 | 17.725 | 0.326 | 2019-2001 |
$AN$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| AN 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 11.624 | -21.346 | 44.594 | 0.464 | 2019-2001 |
$BS
$BS$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BS 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -18.738 | -51.619 | 14.142 | 0.245 | 2019-2001 |
$BS$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BS 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -25.71 | -57.388 | 5.967 | 0.105 | 2019-2001 |
$BS$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BS 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -34.015 | -76.337 | 8.307 | 0.108 | 2019-2001 |
$GR
$GR$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GR 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -48.474 | -89.466 | -7.482 | 0.023 | 2019-2001 |
$GR$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GR 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -6.768 | -59.155 | 45.62 | 0.788 | 2019-2001 |
$GR$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GR 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -19.036 | -54.623 | 16.551 | 0.274 | 2019-2001 |
$pp
$pp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| pp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -64.815 | -100.543 | -29.086 | 0.001 | 2019-2001 |
$pp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| pp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -2.433 | -50.649 | 45.783 | 0.916 | 2019-2001 |
$pp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| pp 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 1.79 | -36.185 | 39.766 | 0.922 | 2019-2001 |
$rf
$rf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| rf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -7.515 | -30.699 | 15.668 | 0.5 | 2019-2001 |
$rf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| rf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -7.886 | -35.016 | 19.245 | 0.545 | 2019-2001 |
$rf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| rf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 10.934 | -15.022 | 36.89 | 0.383 | 2019-2001 |
$wf
$wf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| wf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -41.049 | -84.594 | 2.497 | 0.063 | 2019-2001 |
$wf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| wf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -40.808 | -80.859 | -0.757 | 0.046 | 2019-2001 |
$wf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| wf 2001 vs. 2019 bulk year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -57.661 | -102.561 | -12.761 | 0.015 | 2019-2001 |
NaNs producedNaNs producedNaNs produced
$ANpp
$ANpp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| ANpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -112.068 | -168.378 | -55.758 | 0.013 | 2019-2001 |
$ANpp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| ANpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -41.107 | -126.069 | 43.856 | 0.173 | 2019-2001 |
$ANpp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| ANpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 26.321 | -162.239 | 214.88 | 0.609 | 2019-2001 |
$ANrf
$ANrf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| ANrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -53.528 | -70.509 | -36.547 | 0.005 | 2019-2001 |
$ANrf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| ANrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -20.65 | -71.163 | 29.863 | 0.221 | 2019-2001 |
$ANrf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| ANrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 20.616 | -7.758 | 48.989 | 0.089 | 2019-2001 |
$ANwf
$ANwf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| ANwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -25.061 | -123.937 | 73.815 | 0.39 | 2019-2001 |
$ANwf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| ANwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -18.73 | -25.823 | -11.636 | 0.007 | 2019-2001 |
$ANwf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| ANwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -26.671 | -63.212 | 9.87 | 0.088 | 2019-2001 |
$BSpp
$BSpp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BSpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -41.7 | -75.535 | -7.865 | 0.034 | 2019-2001 |
$BSpp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BSpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -20.628 | -38.25 | -3.006 | 0.037 | 2019-2001 |
$BSpp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BSpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 9.436 | -108.88 | 127.751 | 0.764 | 2019-2001 |
$BSrf
$BSrf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BSrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -69.149 | -82.933 | -55.365 | 0.002 | 2019-2001 |
$BSrf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BSrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -72.152 | -162.05 | 17.746 | 0.075 | 2019-2001 |
$BSrf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BSrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -61.935 | -70.107 | -53.763 | 0 | 2019-2001 |
$BSwf
$BSwf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BSwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -67.026 | -82.623 | -51.43 | 0.002 | 2019-2001 |
$BSwf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BSwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -107.373 | -168.048 | -46.697 | 0.017 | 2019-2001 |
$BSwf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BSwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -139.763 | -241.947 | -37.579 | 0.028 | 2019-2001 |
$GRpp
$GRpp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GRpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -89.459 | -114.521 | -64.397 | 0.004 | 2019-2001 |
$GRpp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GRpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 25.06 | -77.741 | 127.86 | 0.405 | 2019-2001 |
$GRpp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GRpp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 49.63 | -228.753 | 328.013 | 0.524 | 2019-2001 |
$GRrf
$GRrf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GRrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -13.237 | -147.681 | 121.206 | 0.713 | 2019-2001 |
$GRrf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GRrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 69.3 | NaN | NaN | NaN | 2019-2001 |
$GRrf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GRrf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 166.594 | NaN | NaN | NaN | 2019-2001 |
$GRwf
$GRwf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GRwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -55.927 | -101.954 | -9.899 | 0.035 | 2019-2001 |
$GRwf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GRwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -77.344 | -162.104 | 7.416 | 0.059 | 2019-2001 |
$GRwf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GRwf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 0.089 | -36.296 | 36.474 | 0.993 | 2019-2001 |
$AN
$AN$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| AN 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -63.553 | -94.114 | -32.991 | 0.001 | 2019-2001 |
$AN$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| AN 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -26.829 | -48.81 | -4.847 | 0.022 | 2019-2001 |
$AN$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| AN 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 6.755 | -37.035 | 50.545 | 0.738 | 2019-2001 |
$BS
$BS$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| BS 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -59.292 | -79.352 | -39.232 | 0 | 2019-2001 |
$BS$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| BS 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -66.718 | -102.636 | -30.799 | 0.002 | 2019-2001 |
$BS$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| BS 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -64.087 | -123.626 | -4.549 | 0.037 | 2019-2001 |
$GR
$GR$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| GR 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -52.874 | -112.814 | 7.065 | 0.078 | 2019-2001 |
$GR$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| GR 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -7.054 | -105.253 | 91.146 | 0.873 | 2019-2001 |
$GR$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| GR 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 47.18 | -45.261 | 139.621 | 0.278 | 2019-2001 |
$pp
$pp$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| pp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -81.076 | -110.49 | -51.662 | 0 | 2019-2001 |
$pp$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| pp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -12.225 | -41.041 | 16.591 | 0.367 | 2019-2001 |
$pp$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| pp 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 28.462 | -20.711 | 77.636 | 0.226 | 2019-2001 |
$rf
$rf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| rf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -45.305 | -84.505 | -6.104 | 0.028 | 2019-2001 |
$rf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| rf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -23.261 | -119.616 | 73.094 | 0.593 | 2019-2001 |
$rf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| rf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| 9.738 | -101.496 | 120.971 | 0.847 | 2019-2001 |
$wf
$wf$`10`
| 0-10 cm | ||||
|---|---|---|---|---|
| wf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -49.338 | -79.02 | -19.657 | 0.004 | 2019-2001 |
$wf$`20`
| 10-20 cm | ||||
|---|---|---|---|---|
| wf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -67.816 | -112.331 | -23.3 | 0.007 | 2019-2001 |
$wf$`30`
| 20-30 cm | ||||
|---|---|---|---|---|
| wf 2001 vs. 2019 inc year | ||||
| year.diff | year.lwr | year.upr | year.p.adj | Pairs |
| -55.448 | -108.861 | -2.035 | 0.043 | 2019-2001 |
NANA
Call:
lm(formula = d14c_blk ~ Al_ox + lyr_bot + year, data = sra.01.19.min.reps)
Residuals:
Min 1Q Median 3Q Max
-89.532 -27.114 -0.197 26.055 115.938
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 108.5757 6.0801 17.858 < 2e-16 ***
Al_ox -1.9945 0.1687 -11.822 < 2e-16 ***
lyr_bot -3.0592 0.2487 -12.300 < 2e-16 ***
year2019 -22.4127 4.0521 -5.531 6.71e-08 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 36.11 on 314 degrees of freedom
Multiple R-squared: 0.527, Adjusted R-squared: 0.5225
F-statistic: 116.6 on 3 and 314 DF, p-value: < 2.2e-16
The goal of this modeling exercise is to see how parent material and climate/ecosystem affect estimates of soil carbon ages and transit times. Bulk soil 14C observations from 2001, 2009, and 2019 will be used to constrain the carbon models, as well as observations of 14C-CO2 from laboratory soil incubations of soils collected in 2001 and 2019. Previous work has indicated that the carbon stocks at these sites is likely at equilibrium, so we will apply the steady-state assumption to the modeling.
One pool models have been shown repeatedly to be inadequate for describing soil carbon dynamics. However, as simple models are easier to constrain, we will start with a two-pool parallel and two-series models, as these are the simplest model system beyond the single pool approach.
The two-pool parallel model requires the following parameters: * decomposition constants for each pool (k1, k2) * input partitioning coefficient (\(\gamma\)) * steady-state carbon stocks (C) * inputs (I) * initial values of 14C 1 The two-pool series model requires the following parameters: * decomposition constants for each pool (k1, k2) * transfer coefficient (\(\alpha\)) * steady-state carbon stocks (C) * inputs (I) * initial values of 14C
Decomposition rates (k) are related to the amount of 14C in a pre-bomb system (fraction modern, F) at steady-state by the following equations (cf. Schuur, Druffle, and Trumbore, 2016): >Eq. 1
\[F = \frac{k}{k + \lambda}\] >Eq. 2
\[k = \frac{\lambda \cdot F}{1 - F}\] >where \(\lambda\) is the radioactive decay constant (1/8267).
As the decomposition rates will vary, the initial 14C content can be determined dynamically with equation 1.
Carbon stocks are known, while inputs will be estimated and are related to the steady-state conditions by the following equation: >Eq. 3
\[I = (k_{1} \cdot C_{1}) + (k_{2} \cdot C_{2})\] >where C1 and C2 are the carbon stocks of the two model pools.
Both stocks and inputs can be scaled to the known value of the total carbon pool once the steady-state parameters (k1, k2, and \(\gamma\) or \(\alpha\)) have been determined. Pool sizes are a function of the inputs and input partitioning coefficient at steady-state.
A Monte-Carlo Markov chain approach will be used for parameter estimation in combination with an initial optimization algorithm to determine the best set of initial parameters.
Initial model fitting was performed for both model structures using generous parameter ranges [0, 1] for all three parameters (k1, k2, \(\gamma\) or \(\alpha\)). The initial parameter set was found by fitting the models by eye, followed by optimization with the function “modFit” (R package FME), using the Nelder-Mead algorithm. The best set of parameters found by modFit was then used as the input to a Monte Carlo Markov Chain (MCMC), using the function “modMCMC” (R package FME). The number of iterations for the MCMC optimization was set at 5000 intially, with delayed rejection employed to increase efficiency.
The sum of the mean squared error for the best parameter set was slightly lower for the parallel structure than for the series structure. Additionally, the overall mean error of the residuals was also lower for the parallel structure, moderately so for the bulk C observations but substantially so for the respiration observations (in andesite and granite soils in particular).
However, these initial fits yielded unrealistic parameter estimates for multiple sites, particularly at the lower depths. Additionally, the modFit output showed very high correlation between the parameters for both model structures (slightly higher for the two-pool series model).
Optimizing the parameter set requires imposing costs and optionally constraining the allowable range of values for each parameter. Given that we only have data for three time points, this is a relatively sparse data set for constraining these models. Accordingly, the optimization procedure will benefit from a priori constraints of the allowable parameter ranges. For example, since we assume that the system cannot be adequately modeled as a single homogenous reservoir, we will ensure that the optimization procedure cannot collapse the two-pool system into a single pool. This can be mitigated in the two-pool parallel optimization by constraining \(\gamma\) (i.e. the percentage of the inputs entering the fast pool) to a range of 50% to 95%. Similarly, for the two-pool series model structure we can constrain the range of the transfer coefficient to be between 0.0 and 0.1, ensuring that some carbon remains in the fast cycling pool.
Additionally, to enforce a relatively fast cycling pool and relatively slower cycling pool, we will loosely constrain the intrinsic decomposition rates as well (both model structures):
k1: [0.02, 1.00] (50 to 1 year) k2: [0.0001, 0.02] (10,000 to 50 years)
Finally, the models will be run to enforce steady-state, i.e. with unvarying carbon stocks. The amount of carbon observed in the system will be used in the cost function in addition to the radiocarbon observations made in 2001, 2009, and 2019. The inputs will be estimated from net ecosystem exchange (NEE) data measured at nearby eddy covariance sites: Blodgett experimental forest (AmeriFlux), Lower Teakettle (NEON), and Soaproot Saddle (NEON). Alternatively, using correlations between fluxes measured from these eddy covariance towers and GPP estimated from satellite retrievals of SIF, estimates can be made for inputs at the pixels corresponding to each site location.
Error in lapply(mod.fits.2pp3, "[[", 1) :
object 'mod.fits.2pp3' not found